( N AS A-CR- 133888) RELATIVISTIC PARTICLE N73-31092 

MOTION IN NONUNIFORM ELECTROMAGNETIC WAVES 
(California Univ.) 14 p HC S3. 00 

CSCL 20 N Unclas 

G3/07 15439 



UNIVERSITY OF CALIFORNIA 

LOS ANGELES 




Relativistic Particle Motion in Nonuniform 
Electromagnetic Waves 


G. Schmidtt and T. Wilcox 


PPG-152 


June 1973 


Department of Physics 
University of California 

4 

Los Angeles, California 90024 


+ Permanent address: Department of Physics, Stevens Institute of 

Technology, Hoboken, New Jersey 07030. 



Relativistic Particle Motion in Nonuniform 


Electromagnetic Waves 


G. Schmidt^ and T. Wilcox 
Department of Physics 

University of California, Los Angeles, California 90024 


Abstract 

It is shown that a charged particle moving in a strong nonuniform 
electromagnetic wave suffers a net acceleration in the direction of 
the negative intensity gradient of the wave. Electrons will be expelled 
perpendicularly from narrow laser beams and various instabilities can 


result . 
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Particles moving in electromagnetic wave fields strong enough to 
drive them relativistic have been extensively investigated in recent years, 
both with regard to cosmic ray production in pulsars^ and laser particle 
interaction. Linearly polarized plane waves can lead to particle accelera- 
tion in the forward direction while particles in a circularly polarized 
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wave are thought to produce d.c. magnetic fields . Plasma effects 

4 

including instabilities have also been studied in some cases. Some 

problems in which the nonuniformity of the wave plays a role have 
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recently been investigated by Kaw and Kulsrud b and Vittitoe and Wright . 

Here a general treatment is presented of particles moving in weakly 

nonuniform waves where the wave intensity varies slowly on the space and 

time scale of the particle oscillation period. It has long been known that 

nonrelativistic particles suffer a net acceleration in such a wave in the 

7 

direction of the negative intensity gradient . This behavior is characterized 
by a ,l ponderomotive force” and leads to mode coupling and instabilities in 

g 

plasmas . While the lowest order particle motion in the nonrelativistic 
case is that of a harmonic oscillator the relativistic case is more involved, 
but it will be shown that field nonuniformities lead to a generalized 
ponderomotive acceleration tending to expell particles from regions of 
strong field intensity. 

The equation of motion of an electron in an electromagnetic field may 
be written as 
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where is the four velocity and A^ the vector potential . This equation may 
be cast in the more convenient form 
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Consider first in the lowest order the motion of the particle in a uniform 
plane wave. Adopting the Coulomb gauge we have A^ = (aA a , 0, 0); 
a = constant; A^x.i-ct) = A x (n). Eq. (2) reads now in components 
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and from Eqs. (3b) and (3c) 


37 (U„ - cy) = 0 (4) 

One may integrate Eqs. (3a) and (4) to obtain constants of motion. 

Now we let A^ = (aA A , bB(n),0) where a and b vary slowly on scale lengths 
much larger than the particle and wave oscillation periods. From V*^t = 0 

t 3a 9b _ , 9B 

** • * 35C B * b nr = 0 (5) 

The two dominant terms are the first and the last one giving 
bB 'v aA^ where X is the wave length and L the scale length of variation of 
the slowly varying a and b. 

Equation (2) yields now 


4- (1 - — aA.) = 

dx v 1 m 


-It . 3, |L 

m t i 3x. 


37 < u " - <*> ■ -sM* (db'Hx ) * 



-4- 


where terms on the right hand sides represent small perturbations due to 

field nonuniformities and the longitudinal bB terms have been neglected, 

being of order A/L times the perturbation terms arising from aX x . These 

2 22 

two equations plus the exact cy = Sc +U M +U A equation (arising from = 

constant) are the complete set of equations to be solved. 

Substituting now the lowest order solution for iJ from the integral 
Eq. (3a) to the right hand side of Eqs. (6) and (7) and integrating over a 
period of particle oscillation yields 



where A represents the change of a quantity over an oscillation period. This 
period can be characterized by the time it takes for the particle to complete 
an oscillation between equal values of A x . The integral in Eqs. (8) and (9) 
is easily evaluated 

| A^dx = j A 1 2 (dn/dt)" 1 dn = ^ <A ± S (10) 

X 

where dn/dx = U,,-cy = -K is constant to lowest order from Eq. (4), and 
2 

< A 1 > is half the amplitude square for linear polarization and the amplitude 
square for circular polarization. 

Consider first a light pulse broad in the direction perpendicular to its 

propagation such that 3a 2 /3x JL £ 0. Since for such a short pulse f— + - — )a 2 

c 3 t' 

0 all perturbation quantities vanish. As the pulse moves through an originally 
stationary particle, after the pulse passed the particle is left stationary 
(A t x = 0, A0 X = 0, AU„ = 0). 

In the following we will focus our attention to time independent fields 
3a/3t = 0. Since A(cy) = (cy) * (U 1 AU I +U,,AU,,) one may combine Eqs. (8), (9) 
and (10) to find the change of particle energy during an oscillation period 
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l'or K = const (3a/3x tl = 0) one may integrate Eq. (11) over an oscillation 

period to find <Ay> = 0 where < > signifies the average of a quantity over an 

3a 2 

oscillation period. The integration of Eq. (7) shows that K = Kq+0(A ), 

hence corrections in Eq. (11) due to the variation of K give higher order 
terms. Hence one finds in general the interesting result that the average 
particle energy is an adiabatic invariant. 

Introduce now the (proper) time of an oscillation At = j(dn/dx) *dq = 

AK * and use Eqs. (8) and (9) to calculate the particle acceleration over the 


slow time scale. 
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where is the drift velocity of the particle, whose acceleration is due to 
the intensity gradient of the wave. Hence a particle oscillating in a strong 
wave field suffers an acceleration toward the weaker field region and is 
ultimately ejected from the beam. In the process the oscillatory particle 
energy turns into directed energy keeping y constant on the average. Ony may 
multiply Eq. (14) by to find 

2 

2 e 2 2 

Uj + — j <a A a > = constant (15) 

in 


which expresses again the constancy of the sum of oscillatory and drift 

2 2 2 2 

energies. Eq. (15) may be regarded as an energy equation with <f> = e /m <a A ± > 



- 6 - 


to play the role of a potential. A particle when injected from outside the 
beam will be reflected as from a potential barrier. 

Finally we wish to point out some consequences of this acceleration. 
Particles in a laser beam will be accelerated and ejected sideways. The 
presence of a plasma will modify our equations (e.g. the phase velocity of 
waves is no longer c), but the basic effect of ponderomotive acceleration 
is still there. The ejection of particles leads now to a change of dielectric 
function along the beam path and results in self focusing and filamentation 
as in the nonrelativistic case. 

For a circularly polarized wave the lowest order particle motion is 
gyration with = constant, with the electric force providing the 

centripetal acceleration and v*& =0. In the presence of 9a/ 3x the guiding 
center will be accelerated toward the weaker field, as can be seen from Eq. (6) 
with the right hand side providing a constant acceleration. The physical 
reason for this acceleration is easily seen; the electric force acting on 
the particle is stronger along part of its orbit where the intensity is larger 
than on the weaker field side, providing a net accelerating force. The presence 
of a background plasma is known to lead to an axial d.c. magnetic field. This 
may be easily incorporated in the calculation and one finds that the outward 
ponderomotive force leads to an azimuthal drift in a cylindrical beam. A 
perturbation of such an equilibrium can lead to flute type instabilities. 

Details of this problem will be published elsewhere. 
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